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— This  study  investigates  fatigue-creep  interaction 
effects  in  alloys  and  evaluates  the  effectiveness  of 
predictive  models  currently  in  use  by  the  aircraft  engine 
industry.  The  state-of-the-art  crack  growth  rate  prediction 
models  are  supplied  by  the  General  Electric  Company  (MSE) 
model,  and  the  Pratt  and  Whitney  Aircraft  Group  (SINH) 
model.  They  are  used  to  predict  crack  growth  rates  under  a 
range  of  conditions  which  involve  fatigue-creep  interactions. 

Another  aspect  of  this  study  involves  the  development 
of  an  empirical  model  to  predict  fatigue-creep  crack  growth 
based  on  creep  crack  growth  rate  data  and  knowledge  of  the 
loading  wave-form.  This  study  is  primarily  directed  toward 

high  temperature  1000  F)  fatigue-creep  interaction  at  low 

?  - 

test  frequencies  and  positive  stress  ratios. 

The  SINH  model  proves  to  be  more  accurate  than  the  MSE 
model  in  predicting  crack  growth  rates  for  the  data 
analyzed.  Both  models  predict  linear  relationships  for  the 
variation  of  crack  growth  rates  (da/dN)  with  the  length  of 
hold-time  or  the  frequency  rate  on  logrithmic  coordinates. 
The  MSE  model  is  developed  for  AF 115,  whereas  the  SINH  model 
has  the  ability  to  adapt  to  any  high-strength  material. 

The  predictive  model,  developed  in  this  paper,  compares 
well  with  experimental  fatigue-creep  crack  growth  rate  data 
that  are  mostly  in  the  time-dependent  regime,  that  is,  when 
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Evaluation  of  Fatigue-Creep  Crack 
Growth  in  an  Engine  Alloy 


I .  Ininadusiian 

Ever  increasing  performance  requirements  for  U.S.  Air 
Force  turbine  engines  have  placed  stringent  demands  on 
material  capability.  A  need  to  increase  the  ratio  of  engine 
thrust  to  weight  has  forced  design  stresses  to  higher 
levels,  approaching  the  materials  yield  strength.  At  the 
same  time  engine  operating  temperatures  have  also  continued 
to  increase.  In  the  hot  section  of  the  engine,  turbine 
disks  operate  under  severe  thermal  and  mechanical  loading. 
The  variety  of  operating  conditions  within  the  engine 
dictate  diverse  material  capability  requirements.  Thus,  the 
disks  materials  must  exhibit  both  strength  at  high 
temperature  and  creep  resistance. 

In  the  1960’s,  the  primary  design  limitations  for 
turbine  components  were  material  strength  and  creep 
resistance.  However,  by  1975  low  cycle  fatigue  (LCF) 
performance  had  become  the  life  limiting  factor  for  over  75? 
of  the  major  structural  components  in  advanced  engines  such 
as  the  F-100  which  powers  the  Air  Forces  F-15  and  F-16 
aircraft  (1:355).  The  present  design  philosophy  is  based 
entirely  on  the  concepts  of  fatique  crack  initiation  (2,3). 
Any  additional  component  life  due  to  subcritical  crack 
propagation  following  initiation  is  not  considered. 


According  to  the  initiation  philosophy,  1  out  of  1000  of 
these  components  will  have  a  measurable  flaw,  while  the 
remaining  components  have  not  exhausted  their  initiation 
life.  The  difference  between  the  design  life  of  a  disk  and 
the  mean  life-time  of  all  disks  may  be  an  order  of  magnitude 
or  more.  The  present  forced  retirement  policy  requires  the 
elimination  of  999  statistically  sound  disks  out  of  1000. 
Although  this  design  approach  is  extremely  conservative,  it 
is  considered  necessary  for  major  structural  components 
whose  failure  may  have  catastrophic  results. 

However,  turbine  engine  components  are  becoming 
increasingly  costly  to  replace;  moreover,  they  are  made  of 
materials  which  are  of  limited  domestic  supply.  Therefore, 
both  from  an  economic  and  from  a  strategic  viewpoint,  it  is 
desirable  to  fully  utilize  the  available  component  life, 
provided  it  can  be  done  without  sacrificing  safety.  To 
accomplish  this,  an  alternative  approach  to  life  management 
of  engine  components  has  been  proposed.  This  new  philosophy 
bases  component  life  prediction  upon  both  the  initiation  and 
stable  growth  of  cracks  in  LCF  limited  components.  This 
philosophy  recommends  tracking  of  individual  engine 
components  and  that  removal  from  service  decisions  be  based 
on  individual  component  capability  rather  that  on  a  worst 
case  analysis  for  the  entire  population  of  components. 

The  replacement  of  a  safe  life  design  policy  by  one  of 
"fly  to  a  safe  crack"  has  been  termed  retirement  for  cause 
(RFC)  ( 4 ; 5 ) .  This  approach  is  being  considered  for  use  in 
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military  engines  by  the  U.S.  Air  Force.  RFC  is  based  on 
statistical  considerations  that  most  disks  have  no 
detectable  cracks  when  they  reach  their  design  life  and 
would  then  have  some  useful  life  remaining.  At  present, 
each  component  is  retired  from  service  when  it  has  reached 
its  design  life,  whether  or  not  any  cracks  are  found.  Under 
RFC,  components  would  be  kept  in  service  and  inspected  at 
predetermined  intervals  for  fatigue  cracks.  Components 
with  no  detectable  cracks  would  be  kept  in  service  for 
another  inspection  interval.  Conversely,  if  a  crack  were 
detected,  then  that  part  would  be  retired  for  cause. 
Fracture  mechanics  analysis  would  ensure  that  a  crack 
smaller  than  the  detectable  size  would  not  grow  to  a 
critical  size  in  less  time  than  the  specified  inspection 
interval  multiplied  by  some  factor  of  safety,  e.g.  2  or  3. 
This  procedure  provides  several  opportunities  to  detect  a 
crack  before  it  reaches  its  critical  length.  It  is  obvious 
then  that  the  reliability  of  both  the  inspection  techniques 
and  the  crack  growth  predictions  are  vital  aspects  of  RFC. 

Three  primary  issues  must  be  addressed  for  the 
successful  implementation  of  an  RFC  approach  to  fleet 
management.  The  first  one  concerns  the  inspectability  of 
cracks  or  flaws.  A  nondestructive  evaluation  (NDE) 
procedure  which  is  both  reliable  and  accurate  must  be 
available.  The  second  one  deals  with  the  ability  to  predict 
the  rate  of  growth  of  a  crack  under  the  loading  and 
temperature  conditions  which  will  be  encountered  in  the 
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design  mission.  The  third  one  considers  the  economic 
aspects  of  the  RFC  concepts  including  inspection  costs, 
replacement  costs,  and  statistics  of  safety  and  reliability. 
In  this  paper,  the  predicitons  of  crack  growth  behavior  will 
be  considered. 

A  typical  loading  spectrum  for  a  turbine  engine  disk  is 
composed  of  a  mix  of  high  and  low  stress  fatigue  cycles  and 
various  periods  of  sustained  loading.  Typical  fatigue 
cycles  generally  have  a  frequency  of  0.1Hz  or  slower,  and 
are  of  the  order  of  minutes.  At  high  temperatures  crack 
propagation  produced  by  mission  loading  is  quite  complex, 
requiring  the  capability  to  predict  the  effects  of  fatigue 
and  creep,  and  their  interactions.  It  is  the  aim  of  this 
study  to  investigate  these  fatigue-creep  interaction 
effects.  This  is  done  by  implementing  the  state  of  the  art 
crack  growth  rate  prediction  models  currently  used  in 
industry.  These  models  were  developed  by  the  General 
Electric  Company  (6)  and  the  Pratt  and  Whitney  Aircraft 
Group  (7)  to  predict  crack  growth  rate  under  a  range  of 
conditions  which  involve  f a t i gu e- c r ee p  interactions. 
Another  aspect  of  this  study  involved  the  development  of  an 
empirical  model  to  predict  fatigue-creep  crack  growth  based 
on  creep  crack  growth  data  and  knowledge  of  the  loading 
wave-form.  This  effort  is  primarily  directed  toward  high 
temperature  (21000  F)  fatigue-creep  interaction  at  low  test 
frequencies  and  positive  stress  ratios,  as  this  represents 
the  main  area  of  concern. 


The  crack  growth  models  used  were  those  developed  by 
the  General  Electric  Company  (6)  and  the  Pratt  and  Whitney 
Aircraft  Company  (7).  Both  crack  growth  computer  programs 
are  empirical  models  based  on  the  sigmoidal  shape  of  the 
crack  growth  rate  (da/dN)  versus  stress  intensity  range 
(aK)  curve.  An  example  of  such  a  graph  is  illustrated  in 
Figure  1 . 

Both  models  attempt  to  predict  the  crack  growth  rate 
behavior  in  Regions  I  and  III  along  with  the  usual  straight 
line  portion  of  Region  II.  In  many  materials,  especially  at 
high  temperatures  and  stresses,  the  linear  portion  of 
Figure  1  is  very  narrow  or  nonexistent  (8:226).  luerefore, 
accurate  and  reliable  predictions  of  crack  growth  rate 
behavior  in  Regions  I  and  III  are  required  by  these  crack 
growth  models.  The  model  developed  by  the  Pratt  and  Whitney 
Company  uses  a  hyperbolic  sine  equation  fitted  to  a  set  of 
experimental  data  by  the  method  of  least  squares  (SINH 
model).  This  model  has  the  capability  of  predicting  crack 
growth  for  various  temperatures,  stress  ratios,  frequencies, 
and  hold-times  of  various  duration.  The  form  of  the 
hyperbolic  sine  equation  used  in  this  model  is  given  in 
Eq  (1) 


Log  (da/dN)  =  C1  sinh  (C2(logUK)  +  C3))  +  C4 


(1) 


where 
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C-|  =  material  constant  (0.5  for  AF 1 15) 


C2  =  C2(R,  T,  f,  tn) 


c3  =  C3(R,  T,  f,  tH) 
Ci|  =  Cj|(R,  T,  f,  tfl) 


where  Ris  the  stress  ratio,  T  the  temperature,  f  the 
frequency,  and  t^  the  hold-time.  The  model  is  formulated 
in  an  iterative,  interactive  computer  program  for  computa¬ 
tional  purposes.  If  the  program  is  supplied  with  at  least 
two  but  no  more  than  eight  sets  of  data  in  which  only  one 
test  parameter  is  varying  (e.g.  stress  ratio)  then  the  SINH 
model  can  interpolate  for  that  parameter  in  the  range  given 
with  the  data.  For  example,  if  the  program  is  supplied  with 
three  data  sets  having  stress  ratios  of  .1,  .5,  and  .9, 
then  the  SINH  model  will  provide  the  appropriate  constants 
needed  for  interpolating  between  these  stress  ratios.  For 
this  example,  where  all  test  parameters  are  constant  except 
for  stress  ratio  (R),  Cg  through  C4  are  only  functions  of 
stress  ratio.  One  can  then  produce  an  equation  for  a  crack 
growth  rate  curve  as  a  function  of  the  stress  intensisty 
range  aK  for  any  stress  ratio  in  this  range.  This  SINH 
model  will  be  described  in  more  detail  in  the  Numerical 
Analysis  section. 

The  General  Electric  Company  developed  a  model  for 
predicting  crack  growth  based  on  a  modified  sigmoidal 


equation  (MSE).  This  equation  consists  of  the  products  of 
the  exponential  and  natural  logarithmic  functions  raised  to 
constant  powers: 

da/dN  =  exp(B)(AK/ AK*)p(ln(AK/A  K*))Q(ln(AKc/AK))D  (2) 

where 

aK  =  stress  intensity  range 
aK*  =  threshold  aK  from  model 
aKc  =  critical  aK  from  model 

B,P,Q,D  =  constants  for  curve  shape. 

This  model  was  applied  specifically  to  the  material  AF115 
(nickel-base  powder  alloy).  Thus,  this  model  may  not  apply 
as  well  to  other  materials.  The  input  for  using  the  MSE 
model  consists  of  the  various  test  parameters  such  as  the 
stress  ratio,  the  temperature,  the  frequency,  and  the  hold- 
times.  The  constants  for  Eq  (2)  are  the  output. 

Both  models  use  experimental  data  for  AF115  to  evaluate 
their  capability  to  predict  crack  growth  rate  behavior. 
This  material  is  a  nickel-base  powder  super  alloy  developed 
by  the  General  Electric  Company  under  Air  Force  Material 
Laboratory  (AFML)  sponsorship  (9).  A F 1 1 5  is  a  creep 
resistant  high  strength  alloy  with  operating  temperatures  up 
to  1400F.  This  alloy  was  designed  to  be  used  in  high 
temperature  environments  where  space  constraints  are 
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critical.  Applications  that  would  meet  this  criterion  would 
primarily  be  turbine  disks.  AF115  is  then  a  prime  candidate 
for  the  study  of  high  temperature  crack  growth  under  cyclic 
loading  with  interspersed  periods  of  hold-time  loading. 

Equations  (1)  and  (2)  obtained  from  the  supplied  models 
were  numerically  integrated  for  an  estimate  of  the 
specimen's  useful  life.  Results  of  the  integration  were 
represented  as  crack  length  (a)  versus  number  of  cycles  (N). 
The  integration  method  was  a  subroutine  supplied  by  the 
International  Mathematical  and  Statistical  Libraries, 
Inc.  (IMSL)  (10).  The  name  of  the  subroutine  is  DVERK  and 
uses  Runga  Kutta  methods  of  orders  five  and  six.  The 
integrated  data  were  plotted  along  with  experimental  a 
versus  N  data  for  each  loading  condition  considered. 
Conclusions  were  drawn  from  these  plots  to  determine  the 
accuracy  and  reliability  of  each  model. 

The  above  discussion  outlines  the  method  used  to  test 
the  SINH  and  MSE  models  for  crack  growth  accuracy.  Next, 
this  study  examined  each  model's  interpretation  of  crack 
growth  rate  (da/dN)  as  a  function  of  hold-time,  then,  as  a 
function  of  frequency.  Hold-time  variation  will  be  used 
here  for  illustrating  the  procedure.  The  SINH  model  was 
supplied  with  three  data  sets  with  hold-times  of  0,  90,  and 
300  seconds.  It  provided  the  necessary  information  needed 
for  interpolating  between  these  hold-times.  Several  hold- 
times  were  then  selected  so  that  several  crack  growth  rate 
equations  as  a  function  of  the  stress  intensity  range  ( aK) 
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can  be  produced  (see  Eq  (1)).  Fixing  the  value  of  aK 
generated  several  values  of  da/dN.  This  crack  growth  rate 
data  was  plotted  against  the  corresponding  value  of  hold¬ 
time.  This  allowed  examination  of  SINH  interpretation  of 
crack  growth  rate  as  a  function  of  hold-time.  The  same 
procedure  was  used  to  examine  frequency  variation.  The 
method  of  analyzing  the  MSE  model  was  conducted  in  a 
similiar  manner. 

The  MSE  model's  interpretation  of  hold-time  and 
frequency  effects  on  crack  growth  was  examined  in  the 
following  manner.  Again  hold-time  is  used  for  illustration 
purposes.  A  hold-time  was  selected  between  0  and  300 
seconds.  This  hold-time  was  input  into  the  MSE  program 
along  with  the  other  constant  test  parameters.  The 
constants  for  Eq  (2)  are  the  output.  Hence,  the  crack 
growth  rate  equation  for  that  hold-time  is  known.  This 
procedure  was  repeated  for  several  different  hold-times. 
The  resulting  set  of  equations  was  handled  in  the  same 
manner  as  was  done  with  SINH's  equations  described  above. 

Another  area  studied  was  the  ability  to  predict 
fatigue-creep  crack  growth  from  pure  creep  data.  An 
empirical  model  was  developed  that  generated  crack  growth 
rate  (da/dN)  data  from  experimental  creep  crack  growth  rate 
(da/dt)  data  and  knowledge  of  the  cyclic  wave  form.  Two 
wave-forms  were  considered  for  this  model.  Figure  2 
illustrates  these  two  wave-forms.  The  first  is  composed  of 
loading  at  a  constant  rate  up  to  a  specific  stress 


intensity,  holding  that  stress  for  a  period  of  time,  then 
unloading  at  the  same  rate.  The  hold-time  at  constant  load 
was  the  changing  parameter.  The  second  wave-form  consisted 
of  loading  at  a  specified  rate,  then  unloading  at  the  same 
rate.  The  rate  of  loading  and  unloading  was  the  varying 
parameter.  They  were  characterized  using  cycle  time.  Cycle 
time  is  defined  to  be  the  time  required  for  either  wave  form 
to  complete  one  cycle  or  period.  By  varying  the  cycle  time, 
various  values  of  da/dN  were  generated  by  the  model.  This 
da/dN  versus  cycle  time  data  was  then  plotted  along  with 
experimental  data  for  analysis  of  the  model's  validity. 
This  analysis  aided  in  understanding  some  interaction 
effects  between  fatigue  and  creep.  The  data  used  in  this 
model's  development  and  analysis  was  supplied  by  the  Air 
Force  Materials  Laboratory  (AFMLX 1 1 ;  1 2).  This  material  is 
Inconel  718,  a  nickel-base  super  alloy  commonly  used  in 
turbine  engines.  Inconel  718  was  chosen  for  this  model 
development  because  there  was  no  creep  crack  growth  data 
available  from  the  AF115  batch.  The  ability  to  develop  this 
model  and  apply  it  to  life  predictions  is  based  on  the 
following  simplifying  assumptions. 

The  first  assumption  is  that  Linear  Elastic  Fracture 
Mechanics  (LEFM)  correctly  represents  the  stress  and  strain 
conditions  ahead  of  the  crack  tip.  The  equation  for  the 
stress  intensity  factor  range  (aK),  used  for  integrating 
crack  growth  rate  curves,  does  not  include  any  plastic  zone 
corrections.  Lastly,  a  power  law  is  the  assumed  model  for 


the  creep  crack  growth  rate  as  a  function  of  the  stress 
intensity  factor.  In  this  instance,  the  power  law  takes  the 


form  of  Eq  (3) 


da/dt  =  CKm 


where 


empirical  constants  to  be  determined 
from  experimental  data 


K  s  stress  intensity  factor. 


As  presented  in  the  results  and  discussion  section,  this 
assumption  fits  the  data  reasonably  well  for  the  range  of 
stress  intensities  being  considered. 

To  summarize,  two  models  were  used  to  make  specimen 
life  predictions,  conclusions  were  made  on  each  model's 
accuracy  and  consistency  based  on  these  predictions.  The 
models  were  used  to  investigate  their  ability  to  predict  the 
effect  of  hold-time  and  frequency  on  crack  growth  rates. 
Finally,  an  empirical  model  was  developed  to  aid  in 
understanding  the  effect  of  fatigue-creep  interaction. 
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in.  fluiafiridal  Analysis 


Detailed  descriptions  are  presented  here  of  both 
General  Electric's  and  Pratt  and  Whitney's  models.  These 
models  represent  the  state-of-the  art  capability  for  crack 
growth  rate  predictions  for  the  aircraft  engine  industry. 
The  models  are  discussed  in  sufficient  detail  in  order  to 
present  each  model's  theoretical  foundation  and  basis  for 
usage.  Detailed  development  for  these  models  are  given  in 
(6;7).  This  section  also  presents  the  development  of  a 
predictive  model  to  evaluate  the  ability  to  use  creep  data, 
along  with  knowledge  of  the  wave-form,  to  predict  fatigue- 
creep  crack  growth  rates. 
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The  interpolative  hyperbolic  sine  model  (SINH)  is  a 
mathematical  tool  which  provides  crack  propagation  rate  as  a 
function  of  various  parameters.  Pratt  and  Whitney  selected 
the  hyperbolic  sine  equation  to  represent  crack  growth  rate 
(da/dN)  as  a  function  of  the  stress  intensity  range  (aK)  for 
the  following  reasons: 

1.  It  exhibits  the  overall  sigmoidal  shape  of  typical 
da/dN  versus  aK  plots  as  shown  in  Figure  1 . 

2.  All  or  part  of  the  equation  may  be  used  to  fit 
data.  Also,  the  slope  at  the  inflection  can  be 
adjusted  with  the  fitting  constants. 

3.  The  hyperbolic  sine  is  not  periodic  nor 
asymptotic  (e.g.  trigonometric  tangent);  therefore 
when  extrapolation  becomes  necessary,  the  sinh 
behaves  well  at  distances  removed  from  data. 


^  T  « 


V 


r 


The  SINH  model  is  based  on  the  hyperbolic  sine  function  and 
is  given  as 

log(da/dN)  =  ClSinh(C2(log( aK)  +  C3))  +  C4  (4) 

where  the  constants  have  been  shown  (15;  16;  17)  to  be 
functions  of  the  test  frequency,  the  stress  ratio,  and  the 
temperature  for  pure  cyclic  loading: 


C-j  =  material  constant  (0.5  for  AF115) 
C2  =  C2(f»R»T)  =  shape  factor 


C o  =  C-j(f,R,T)  =  abscissa  of  the  inflection  point 
0  J  (-logAK) 

c4  =  C4 ( f , R,  T)  =  ordinate  of  the  inflection  point 

(log(da/dN) ) 


where  f  is  the  frequency,  R  the  stress  ratio  and  T  the 
temperature.  For  interpolating  in  regions  where  there  are 
no  data,  the  constants  C2-C4  are  assumed  linear  functions  of 
the  crack  growth  parameters  as  specified  below 


C2  =  a  +  beFNCTN 


Cn  =  c  +  d#FNCTN 


C3  =  e  +  f#C4 

FNCTN  is  a  function  of  the  individual  test  parameter  which 
is  being  regressed.  Table  1  defines  these  relationships. 
These  relationships  allow  development  of  an  interpolative 
model  of  crack  propagation  as  a  function  of  any  one  test 
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parameter.  Because  of  the  simple  relationships  observed 
between  the  coefficients  of  the  SINH  model  and  the  basic 
crack  propagation  controlling  parameters,  interpolations  are 
straight  forward.  It  is  here  that  the  usefulness  of  the 
model  is  best  demonstrated;  the  SINH  model  provides 
descriptions  of  crack  propagation  characteristics  even  in 
regions  where  data  are  unavailable.  Appendix  A  gives  a 
basic  development  of  the  SINH  model.  For  a  detailed 
development  see  (7). 


Ike  General  fileairisls  UJSEI  ftedel 

The  development  of  General  Electric’s  model  is  based  on 
the  introduction  of  a  sigmoidal  equation  which  had  the 
flexibility  of  correlating  the  complete  range  of 
conventional  cyclic  crack  growth  rate  versus  stress 
intensity  range  data.  The  interpolative  model  was  developed 
using  a  modified  form  of  a  sigmoidal  equation  (MSE)  given  in 
Eq  (2).  The  lower  and  upper  asymptotes  are  expressed  by  aK* 
and  aKc,  respectively.  The  coefficients  Q  and  D  are  shaping 
coefficients  that  control  the  lower  and  upper  section  of  the 
sigmoidal  curve  (see  Figure  3).  Decreasing  the  absolute 
value  of  D  or  Q  results  in  a  sharper  transition  at  the 
appropriate  end  of  the  curve.  The  coefficient  P  in  the  MSE 
allows  control  of  rotation  at  the  inflection  point  of  the 
crack  growth  curve.  Changing  the  B  coefficient  results  in  a 
vertical  movement  of  the  curve.  The  horizontal  location  of 
the  inflection  point  is  determined  by  a  mathematical 
combination  of  aK*,  AKc,  D,  and  Q  (6:37). 
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The  six  coefficients  (B,  P,Q,D,  aK*,  and  aKc)  are  related 
to  temperature,  stress  ratio,  frequency  and  hold-time.  Some 
of  the  coefficients  interact  with  only  one  or  two  of  the 
four  test  variables.  Others,  such  as  the  ones  used  for 
locating  the  inflection  point  are  more  complex  and  are 
related  to  all  four  test  variables  (6:37).  In  many 
situations  during  the  development  of  this  model  by  General 
Electric,  the  coefficients  governing  the  crack  growth  curve 
were  estimated  since  tests  were  not  conducted  at  all 
possible  combinations  of  the  four  test  variables.  The 
estimates  were  achieved  by  considering  the  surrounding  test 
results  in  which  only  one  test  variable  was  different  from 
that  being  examined.  Details  of  each  of  the  derived 
relationships  are  given  in  (6:38-51). 

The  derivation  of  the  model  was  based  essentially  on 
examining  how  the  da/dN  versus  aK  curves  changed  with  a 
change  in  test  parameter.  The  experimental  data  were  then 
incorporated  into  the  MSE  model  through  the  six 
coefficients.  These  coefficients  are  functions  of  the  test 
parameters;  therefore,  changing  the  test  parameters  would 
change  the  six  constants  which  in  turn  would  alter  the  shape 
of  the  MSE.  Modeling  of  the  upper  and  lower  asymtotes  will 
be  presented  as  examples  of  the  derivation  technique.  These 
derivations  are  taken  from  (6:38-39)  and  are  based  on  AF115. 

The  final  crack  length  was  measured  from  each  of  the 
failed  specimens  and  the  stress  intensity  range  was 
calculated  based  on  the  crack  length  and  the  applied  loads. 


No  correlation  of  the  final  stress  intensity  range  to 
frequency,  hold-time,  or  temperature  was  found.  On  the 
other  hand,  the  final  stress  intensity  range  did  depend  on 
stress  ratio;  this  relationship  may  be  written  as  follows: 

aKc  =  AKmax(l-R)  (6) 

The  value  of  AKmax  was  calculated  to  be  122  ksi  (i ri)1  / 2  . 

The  lower  asymptote  aK*,  was  estimated  from  each  set  of 
da/dN  versus  aK  data.  It  was  found  to  be  dependent  on  the 
stress  ratio  and  the  temperature  and  independent  of  the 
frequency  and  the  hold-times.  For  the  stress  ratio  of  0.1 
and  temperatures  between  1000  and  1400  F,  K*  was  estimated 
as  10  ksi  @.n)1/2  .  When  plotted  versus  (1-R)  on  logarithmic 
coordinates,  as  done  by  General  Electric  (6:39),  linear 

relationships  were  observed  for  each  temperature 

(see  Figure  4)  so  that: 

aK*  =  AK*max(1-R)2  (7) 

For  each  temperature,  a  K*max  is  the  intercept  at  the  stress 
ratio  of  zero,  and  n  is  the  slope.  They  were  related  to 
temperature  by: 

n  =  +  b-|  ( T  -  1000)d 

* 

AKmax  =  a2  +  b2(T-1000)e 

where  ai,  a2,  b-|,  b2»  d,  and  e  are  constants  determined  by 
simple  regression  analysis  (Figure  4). 


A  summary  of  the  coefficients  influenced  by  the  four 
test  parameters  are  given  by  General  Electric  (6:53)  and 
reproduced  here  in  Table  2.  These  coefficients  are  used  to 
evaluate  the  constants  in  the  MSE  model.  AKi  and  da/dNi  are 

A 

used  to  evaluate  B.  Also,  da/dNi  was  used  to  find  the  slope 
of  the  curve  at  the  inflection  point  (P)  (6:37).  The  upper 
shaping  parameter  (D)  is  not  listed  in  Table  2,  because  it 
is  calculated  knowing  Q,  AKi,  aKc1  and  aK*  (6:37).  Complete 
derivations  of  these  six  coefficients  are  given  in  (6). 
Appendix  B  presents  the  computer  code  for  the  MSE  model  used 
in  this  study.  This  code  is  an  updated  version  of  the  code 
listed  in  the  appendix  of  (6). 


Zredistiye 

Presented  here  is  the  development  of  a  predictive  model 
to  examine  the  ability  to  predict  fatigue-creep  crack  growth 
from  creep  data  and  knowledge  of  the  wave-form.  Modeling  of 
the  experimental  creep  data  was  the  first  step  of 
development.  The  creep  crack  growth  rate  data  used  in  this 
analysis  is  shown  in  Figure  5.  For  the  stress  intensity 
ranges  used  in  this  study,  creep  crack  growth  rate  data  from 
Inconel  718  can  be  represented  by  a  linear  relationship  on 
logrithmic  coordinates.  The  equation  which  expresses  creep 
crack  growth  rate  as  a  function  of  stress  intensity  factor 


da/dt  =  CKra 
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Figure  5. 

Maximum 

Stress  Intensity  Factor  Kmax, 

versus 

the  Creep  Crack  Growth  Rate  da/dt. 

for  Inconel  718  at  1200  F. 

where 


K  =  stress  intensity  factor 
m  =  slope  of  fitted  line 
C  =  conversion  constant 


The  constants  C  and  m  are  determined  by  fitting  a  straight 
line  to  the  baseline  data  in  Figure  5  over  the  range  of  K 
values  being  considered. 

Two  loading  wave-forms  were  selected  to  evaluate  the 
analytical  model  and  are  presented  in  Figure  6.  The 
frequency  of  the  ramps  and  the  magnitude  of  the  hold-times 
are  variable.  Wave-form  (b)  is  equivalent  to  wave-form  (a) 
when  the  hold-time  is  zero,  and  can  be  considered  to  be  a 
subset  of  (a).  The  first  wave-form  is  divided  into  three 
parts  as  shown  in  Figure  6.  This  wave-form  may  be  expressed 
as  a  function  of  time  and  the  loading  parameters.  These 
equations  are  given  as  follows 


for  part  1 : 


K(t)  =  AKr/( 1  —  R )  +  2AKf^t  from  t=0  to 

t=1/2fx  (10) 


for  part  2: 


K(t)  =  aK/(1-R)  =  constant  (11) 

for  part  3: 

K(t)  =  aK/(1-R)  -  2 aKf ut  from  t=0  to 

t= 1/2  fu 


(12) 


where 


=  Kmax  -  Krain  =  stress  intensity  range 
=  Kmin/Kmax  =  stress  ratio 

=  loading  rate  (cycles/second) 

=  unloading  rate  (cycles/second) 

=  hold-time  (seconds) 


Eqs  (10)  —  (12)  are  substituted  into  Eq  (9)  and  analytically 
integrated  for  one  complete  cycle,  namely  for  part  1: 


^/2fl 


C(aKR/(1-R)  +  2AKf1t)mdt 


for  part  2: 


C( aK/( 1-R) )mdt 


for  part  3: 


■I. 


1  /  2  f , 


C( AK/( 1-R)-2  AKf ut)mdt 


If  m  i  -  1  then  integration  yields 


for  part  1 : 


ai  =  C 


UKR/CI-R)  +  AK)»"+1  ..  (AKR/ ( 1  -R)_)_m+ 1 
2(m+1)AKf1 


for  part  2: 


a2  =  CUK/(1-R))m  tH  (17) 

for  part  3: 

pTaK  -  a  K/(  1-R)  )m+1  +  (AK/(1-R))ra+il  ,iQ, 

a2  =  C -  (1o) 

I  2(m+1  )a  Kf  u 

The  result  of  the  integration  is  the  amount  of  crack  growth 
in  one  loading  cycle.  This  is  defined  as  the  fatigue-creep 
cyclic  crack  growth  rate 

da/dN  =  a1  +  a2  +  a3  (19) 

Equation  (19)  can  be  used  to  predict  fatigue-creep  crack 
growth  rates  from  creep  data  and  knowledge  of  the  wave-form. 
This  equation  was  programmed  into  a  computer  whose  inputs 
are  the  creep  data,  and  the  loading  wave-form  parameters 
defined  in  Eq  (12).  The  program  fits  a  straight  line  to  the 
creep  crack  growth  data  to  produce  the  constants  C  and  m  . 
With  these  constants  and  the  wave-form  parameters,  the 
program  will  return  a  value  for  the  cyclic  crack  growth  rate 
(da/dN).  Calculations  were  conductd  for  a  range  of  values 
of  the  wave-form  parameters.  The  calculated  values  of  da/dN 
were  compared  to  experimental  data  on  Inconel  718.  The 
comparisons  are  presented  in  the  Results  and  Discussion 


Section. 


IV.  Results  and  Discussion 

Data  were  collected  on  AF115  from  a  contract  sponsored 
by  the  Air  Force  Wright  Aeronautical  Laboratories  (6).  All 
cyclic  crack  growth  experiments  were  conducted  at  General 
Electric,  Aircraft  Engine  Group  (AEG),  in  two  servo- 
controlled  electrohydraulic  test  systems  operated  under 
closed-loop  control.  Details  on  specimen  fabrication, 
precracking,  heating,  crack  measuring,  and  data  reduction 
are  given  in  Appendix  C.  For  more  information  on  this 
material  see  (6:2-23). 

Results  of  crack  growth  comparisons  of  the  MSE  and  SINH 
models  are  shown  in  Figures  7  through  9.  Each  graph  is 
composed  of  integrated  crack  growth  rate  curves  from  each 
model  and  experimental  data  for  that  loading  condition. 
Notice  that  for  clarity  the  abscissa  of  Figure  7  is  expanded 
with  respect  to  Figures  8  and  9.  The  crack  growth  rate 
curves  along  with  the  experimental  data  used  for  curve 
fitting  are  plotted  in  Figures  10  and  11.  The  equations 
used  to  produce  these  curves  are  given  below  each  graph. 
The  equations  for  the  SINH  model  shown  in  Figure  10  were 
arrived  at  by  producing  an  interpolative  model  for  hold¬ 
time.  Hold-times  of  0,  90,  and  300  seconds  were  chosen  and 
the  constants  C2  through  C4  were  calculated.  The  process 
employed  by  the  SINH  model  to  arrive  at  this  interpolative 
model  was  discussed  in  the  Numerical  Analysis  Section.  The 
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ison  plot  of  AF115  for  R  =  .1,  T  =  1200  F,  f  =  .25  Hz,  and  tu  =  0  sec. 


Experimental  Results 


Comparison  plot  of  AF115  for  R  =  .1,  T  =  1200  F,  f  =  .25  Hz,  and  t  -  90  sec. 
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Figure  jj.  Stress  Intensity  Range  AK,  versus  Crack  Growth  Rate  da/dN 
using  MSE  Model  for  AF115  Data.  R  =  .1,  T  =  1200  F,  and 


constants  used  in  the  MSEs  model  presented  in  Figure  11  are 
obtained  by  substituting  the  test  parameters  into  the  MSE 
computer  program.  The  test  parameters  used  are  presented 
with  each  figure.  Note  that  as  the  hold-time  increases  from 
0  seconds  (Figure  7),  to  90  seconds  (Figure  8),  and  to  300 
seconds  (Figure  9)  the  number  of  cycles  to  failure 
decreases;  also  note  that  the  predicted  curves  begin  to 
deviate  from  the  experimental  data.  Table  3  illustrates 
this  deviation  by  presenting  the  ratio  of  the  predicted 
number  of  cycles  to  failure  (P)  to  the  actual  number  of 
cycles  to  failure  (A).  A  value  for  P/A  less  than  1.0  indi¬ 
cates  a  conservative  estimate  whereas  a  value  greater  than 
1.0  suggests  an  unconservative  estimate  of  crack  growth. 

Besides  investigating  each  model's  capability  to  pre¬ 
dict  crack  growth  curves,  each  model  was  also  examined  for 
its  ability  to  interpret  the  variation  of  the  crack  growth 
rate  with  a  changing  test  parameter.  Two  test  parameters 
were  considered;  hold-time  and  frequency.  Considering  hold¬ 
time  variation  first,  Figures  12  and  13  represent  each 
model’s  interpretation  of  crack  growth  rate  variation  with 
hold-time.  The  stress  intensity  factor  range,  j  K,  along 
with  the  other  test  parameters  and  values  of  da/dN  were 
computed  for  various  hold-times  from  each  model.  The  method 
used  to  compute  this  data  is  presented  in  the  Approach 
section.  The  test  parameters  held  fixed  are  the  stress 
ratio,  frequency,  and  temperature  with  values  of  0.1, 
0.25  Hertz,  and  1200  F  respectively.  The  two  aK  values 


Table  3.  Comparison  of  Predicted  to  Actual 
Number  of  Cycles  to  Failure. 


Model 

Hold-time 

Predicted  N 

( seconds) 

Actual  N 

P.  &  W. 

0.0 

1.1093 

90.0 

1 . 0405 

300.0 

0.7913 

G.  E. 

0.0 

1.0855 

90.0 

0.814-6 

300.0 

0.5338 

Growth 
loads . 


selected  were  22.77  and  32.76  ksi  (in)’'/2.  The  wave-form 
described  by  these  test  conditions  is  illustrated  in 
Figure  6a.  The  SINH  model  shows  some  nonlinearity  for  the 
lower  load  (aK)  (Figure  12)  whereas  the  MSE  model  appears 
linear  for  both  loads  (Figure  13).  Considering  now  the 
frequency  variation  and  using  the  same  procedure  to  compute 
data  as  with  the  hold-time  variation,  the  results  shown  in 
Figures  14  and  15  were  obtained.  A  hold-time  of  0  seconds 
was  chosen  in  this  case.  All  other  test  parameters  are 
equivalent  to  the  hold-time  test.  The  wave-form  described 
by  these  test  conditions  is  presented  in  Figure  6b.  In 
Figure  14,  the  SINH  model  presents  a  very  nonlinear  inter¬ 
pretation  of  crack  growth  rate  variation  with  frequency  for 
the  data  used.  This  model  seems  to  suggest  that  at  a  aK  of 
22.77  Ksi  (in)-'/2  the  crack  growth  rates  actually  decrease 
with  some  values  of  frequency.  The  MSE  model's 
interpretation,  illustrated  in  Figure  15,  consists  of 
parallel  horizontal  lines  for  the  two  loads  considered. 
Figures  16  and  17  present  experimental  data  on  AF115  and  the 
curves  fitted  to  this  data  by  each  model.  AF115  data  with 
frequency  varying  between  each  data  set  was  incorporated 
into  each  model  so  that  interpolative  results  on  frequency 
variation  could  be  obtained. 

For  comparison  with  actual  data  on  another  material, 
Figure  18  shows  experimental  data  on  Inconel  718.  This 
figure  was  produced  at  the  Air  Force  Materials 
Laboratory  (11).  This  material  was  tested  under  the  same 
conditions  as  AF115  whose  computer  generated  data  are 
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Figure  17*  Stress  Intensity  Range AK,  versus  Crack  Growth 
Rate  da/dn,  using  the  MSE  Model  for  AF115  Data. 
R  =  .1,  T  -  1200  F,  t  =  0.0  sec,  f  =  2. 5-. 025 


represented  by  Figures  12  through  15.  Figure  18  illustrates 
the  experimental  crack  growth  rate  variation  with  cycle  time 
for  nickel-base  super  alloy  718.  The  dashed  and  solid 
curves  in  Figure  18  correspond  to  the  loading  wave-forms 
shown  in  Figure  6a  and  6b  respectively.  This  experimental 
data  was  used  to  examine  the  prediction  abiity  of  Eq  (19). 
The  results  of  this  analysis  follow. 

Figure  19  illustrates  the  accuracy  of  Eq  (19)  in 
predicting  fatigue-creep  crack  growth  rates.  Figure  5 
presents  the  creep  data  on  Inconel  718  used  to  determine  the 
constants  C  and  m  in  Eq  (19).  The  baseline  data  were 
selected  for  curve  fitting.  Figure  19a  presents  the  results 
of  using  Eqs  (19)  and  the  loading  wave-form  shown  in 
Figure  6a.  Figure  19b  presents  the  results  of  using  Eq  (19) 
and  the  loading  wave-form  shown  in  Figure  6b.  For 
Figure  19b,  Eq  (19)  and  experimental  data  are  almost 
identical  for  cycle  times  greater  than  5  seconds.  In 
Figure  19b  the  upper  straight  line  was  obtained  by 
integrating  the  entire  associated  wave-form.  The  lower 
straight  line  is  the  result  of  only  integrating  half  the 
wave— form  shown  in  Figure  19b.  Note  that  the  lower 
straight  line  is  a  better  approximation  to  the  experimental 
data  shown  in  Figure  19b.  Also  note  that  the  break  in  the 
bilinear  lines  from  Figure  19a  to  19b  is  shifted  to  the 
right  by  about  an  order  of  magnitude.  This  change  in  slope 
corresponds  to  a  transition  from  predominantly  fatigue  crack 
growth  to  mostly  creep  crack  growth. 


Total  Cycle  Time,  s 


Figure  18.  Fatigue-Creep  Crack  Growth  Rate  da/dN,  for  the 
two  wave- forms  shown. 


da/dN  (in/cycle)  da/dN  (in/cycle) 


V.  Conclusions  and  Recommendations 


Based  on  the  analytical  and  experimental  results,  the 
following  conclusions  can  be  made  for  the  evaluation  of 
fatigue-creep  crack  growth  of  an  engine  alloy  and  the  models 
used  to  analyze  their  behavior. 

SIM 

The  SINH  model  is  a  flexible  program  that  can  be  used 
to  accurately  predict  crack  growth  rates  for  more  than  one 
material.  All  or  part  of  the  model's  curve  can  be  used  by 

properly  selecting  the  equation's  constants,  C2  through  Ci|* 
Since  the  SINH  curve  does  not  have  asymptotes,  this  equation 
behaves  well  for  extrapolating  to  stress  intensity  ranges 
(aK)  outside  of  those  given  in  the  experimental  data. 

As  shown  in  Figure  14,  SINH's  description  of  crack 
growth  variation  with  increasing  hold-times  was  nonlinear  on 
logarithmic  coordinates.  This  nonlinearity  was  believed  to 
be  caused  by  the  response  of  A F 1 1 5  to  the  change  in  hold¬ 
time.  For  instance,  in  Figure  12,  there  was  a  smaller  shift 
between  the  90  and  300  second  data  sets  than  between  the  0 
and  90  second  data  sets.  This  weak  dependency  on  hold-time 
at  the  longer  duration  cycletime  shows  up  in  Figure  14  as  a 
nonlinearity  for  cycle  times  greater  than  90  seconds.  At  a 
load  of  22.77  ksi  in  1/2.  For  hold-time  variation,  the  SINH 
model  assumes  that  the  constants  C2-C4  in  Eq  (4)  vary 
according  to  Eq  (5),  where  FNCTN  is  equal  to  -log(tH).  The 
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method  of  determining  the  interpolative  constants  a  through 
f  in  Eq  (5)  is  described  in  Appendix  C.  If  the  experimental 
data,  used  to  determine  a  through  f,  does  not  correlate  with 
this  assumed  model  then  the  model  will  produce  nonlinear 
results,  as  shown  in  Figure  14.  Since  there  is  a  weak 
dependency  on  hold-time  for  hold-times  greater  than 
90  seconds  (Figure  12),  the  SINH  model  produced  nonlinear 
results  in  Figure  14  for  cycle  times  greater  than 
90  seconds.  Had  the  experimental  data  correlated  with  the 
assumed  model  in  SINH,  the  results  of  Figure  14  should  be 
similar  to  Figure  15.  Notice  that  in  Figure  14  at  a  a K  of 
32.76  ksi  (in)1/2  that  the  result  is  more  linear.  Referring 
to  Figure  12,  there  was  a  wider  shift  between  the  90  and  300 
second  data  sets  at  a  aK  of  33  ksi  ( in )1/2.  This  confirms 
that  as  the  dependency  on  hold-time  strengthens  the  SINH 
model  will  produce  more  linear  behavior. 

All  data,  in  Figure  18,  can  be  fitted  to  a  single  SINH 
curve.  Therefore,  the  crack  growth  rate  (da/dN)  do  not  vary 
with  frequency.  An  interpolative  model  for  frequency  was 
developed  from  the  SINH  model  to  examine  SINH's  description 
of  crack  growth  rate  with  changing  frequency.  This  inter¬ 
polative  model  should  produce  no  variation  of  crack  growth 
rates  with  changing  frequency.  But  as  seen  in  Figure  16 
this  was  not  the  case.  In  Figure  18,  the  inflection  points 
(indicated  by  the  solid  squares)  for  each  of  the  three  SINH 
curves  shift  up  and  to  the  right  as  the  frequency  decreases. 
The  three  data  groups  that  these  curves  were  fitted  to  also 


shift  in  the  same  direction.  Because  of  this  shift  in  the 
experimental  data,  the  SINH  model  does  not  produce  the  same 
equation  for  the  three  data  groups.  Hence  Figure  16  is  not 
independent  of  frequency.  The  same  argument  can  be  made  for 
the  nonlinearity  in  Figure  16  as  was  made  for  Figure  1M  in 
the  previous  paragraph. 

The  MSE  model  only  requires  the  four  test  parameters 
(stress  ratio,  temperature,  frequency,  and  hold-time)  to 
arrive  at  the  associated  crack  growth  rate  equation.  Again, 
this  model  is  only  applicable  to  AF115  material.  This  model 
has  the  ability  to  fit  experimental  data  which  is 
nonsymmetric.  In  other  words,  the  model  is  not  required  to 
be  symmetric  about  its  inflection  point.  Models  regarding 
the  interpretation  of  crack  growth  rates  with  hold-time  or 
frequency  variation  were  incorporated  into  the  MSE  model. 
These  models  were  intended  to  be  linear  in  logarithmic 
coordinates  (see  Figures  13  and  15). 

£IHtt=143E  £fiinj2acisen 

The  examination  of  the  SINH  and  the  MSE  models'  ability 
to  predict  crack  growth  rates  has  led  to  the  following 
conclusions.  First,  both  models  lose  accuracy  as  hold-time 
increases.  Second,  the  number  of  cycles  to  failure  decrease 
with  increasing  hold-times.  Lastly,  on  the  average  the  SINH 
model  produced  more  accurate  P/A  values  than  the  MSE  model 
for  the  data  analyzed  (see  Table  3). 
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The  validity  of  Eq  (19)  was  strictly  within  the  creep 
dominant  range  of  the  experimental  data  used  for  comparison. 
This  model  produced  accurate  predictions  when  the  frequency 
was  greater  than  5  seconds  or  the  frequency  was  less  than 
.02  Hertz.  Therefore,  this  model  should  only  be  implemented 
within  these  regions  when  investigating  the  crack  growth 
rate  behavior  of  Inconel  718.  From  Figure  19,  even  for  very 
low  frequencies  (<  .002  Hz),  most  of  the  damage  accrued 
during  cycling  occurs  only  over  50%  of  the  wave-form.  This 
portion  of  the  wave-form,  over  which  damage  occurs,  may  be 
the  loading  portion,  unloading  portion,  or  a  mixture  of 
both. 

A  sufficient  body  of  AF115  data  should  be  collected  in 
order  to  produce  a  graph  similar  to  the  one  shown  in 
Figure  18.  This  new  information  on  AF115  can  be  used  to 
determine  if  either  SINH’s  or  MSE’s  interpretation  of  crack 
growth  rate  with  cycle  time  is  correct  for  the  data  analyzed. 

As  flexible  as  the  SINH  model  is,  engineering  judgement 
must  be  called  upon  when  analyzing  data.  Seldom  will  the 
model  produce  the  best  possible  results  on  the  first 
attempt.  One  must  examine  the  data  beforehand  to  determine 
the  appropriate  input  parameters  required  to  arrive  at  the 
optimum  output.  The  procedure  used  to  arrive  at  this  best 
fit  to  a  set  of  data  is  more  likely  to  be  iterative,  where 


previous  tries  are  revised  in  hopes  of  obtaining  better 
curve  fits  to  experimental  data.  The  SINH  model  is  a 
sophisticated  curve  fitting  and  data  regression  program  that 
is  specialized  for  cyclic  crack  growth  rate  data.  More  than 
a  few  weeks  are  required  to  fully  understand  all  aspects  of 
this  model. 

Enhanced  understanding  of  fatigue-creep  processes  would 
be  achieved  if  an  empirical  or  analytical  model  could  be 
developed  to  explain  the  crack  growth  behavior  at  cycle 
times  less  than  50  seconds  for  the  experimental  data  given 
in  Figure  19b.  This  region  consists  of  mixed  forms  of 
fatigue  and  creep  crack  growth  processes. 
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The  hyperbolic  sine  is  defined  as 

y  =  sinh  x  =  (exp(x)  -  exp(-x))/2  (A.1) 

and  when  presented  in  cartesian  coordinates  it  appears  as 
shown  in  Figure  20.  The  function  is  zero  at  x=0  and  has  its 
inflection  there. 

The  introduction  of  the  four  regression  coefficients 
Ci  through  C4  permits  relocation  of  the  point  of  inflection 
and  scaling  of  both  axes.  In  the  equation 

(y-C4)  =  SINH  (x+C3)  (A. 2) 

C3  establishes  the  horizontal  location  of  the  hyperbolic 
sine  point  of  inflection,  and  C4  locates  its  vertical 
position. 

To  scale  the  axes,  and  C2  are  introduced 

(y-C4)/C-|  r  sinh  (C2(x+C3))  (A. 3) 

which  can  be  written  as 

y  =  C,  sinh  ( C2( X+C3 ) )  +  C4  (A. 4) 

of  which  Eq  (4)  is  a  special  case  where  y  =  log  (da/dN)  and 
x  =  log  (aK).  Note  that  and  C2  are  dimensionless  and  can 
be  thought  of  as  stretching  the  curve  vertically  and 
horizontally  respectively.  Test  results  have  shown  that, 


for  a  given  material,  Ci  can  be  fixed  without  adversely 
affecting  model  flexibility  (7:36-100). 
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Interpolative  modeling  of  crack  propagation  as  a 
function  of  operating  parameters  such  as  holdtime,  frequency, 
stress  ratio,  and  temperature  requires  a  multiple  regression 
capability  which  allows  simultaneous  consideration  of 
several  different  data  sets,  each  differing  from  the  others 
by  only  one  operating  parameter.  The  model  is  to  have  broad 
interpolative  capability  with  behavior  at  one  condition  used 
to  describe  crack  propagation  at  another  condition,  since 
all  possible  stress  conditions  can  not  be  tested  in  the 
laboratory . 

Pratt  and  Whitney  Aircraft  has  developed  a  mathematical 
technique  to  accomplish  this.  Individual  sets  of  data  are 
treated  independently  relative  to  some  of  the  SINH 
coefficients,  while  the  entire  collection  is  treated  as  an 
entity  with  respect  to  the  interpolative  coefficients 
(C33-C38  discussed  later). 

The  method  of  least  squares  is  the  computational 
procedure  employed  in  the  modeling  of  Eq  (A.M)  to  the  data. 
The  goal  of  this  procedure  is  to  determine  model 
coefficients  for  which  the  resulting  curve  through  the  data 
will  have  the  minimum  summed  squared  error  between 
calculated  and  observed  values  for  the  dependent  variable 
(see  Figure  (21)).  In  this  instance,  the  independent  and 
dependent  variables  are  log  ( aK)  and  log  (da/dN), 
respectively . 


x  =  Log(AK) 


Error  Determination  between 
Observed  and  Calculated  Data 


The  sum  of  the  squared  errors  is  defined  as 


where  N  is  the  number  of  data  points. 

Since  YCALi  =  f(C2,  C3,  C4,  X*)  from  Eqs  (A. 4) 
function  of  C2>  C3,  and  C4. 

Now,  E2  will  be  a  minimum  when  each  of 
derivatives  is  zero  simultaneously.  That  is 

sE  2 
3  C2 

3E2 

7C3 

3_e£ 

3  C4 

when  f  is  the  SINH  model, 

Ei  =  C<|  sinh  ( C2< X±  +  C3 ) )  +  C4  -  Yi 

and 


2E3  E 
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=  C1  cosh  (C2(Xi  +  C3))(Xi  +  c3) 

3  E 


(A. 5) 


E  is  also  a 

its  partial 


(A. 6) 


(A. 7) 


(A. 8) 
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(A. 10) 
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Ct  cosh  (C2(Xi  +  C3))(C2) 


(A.1 1 ) 
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(A. 1 2) 
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Now,  substituting  Eqs  (A. 9)  -  (A. 12)  into  Eqs  (A. 6),  (A. 7), 
and  (A. 8),  and  solving  the  resulting  three  simultaneous 
equations  provides  values  for  C2,  C3,  and  C4  for  which 
Eq  (A. 5)  will  be  a  minimum.  In  this  instance,  the  resulting 
simultaneous  equation  are  nonlinear  in  C2>  C3,  and  C4.  The 
solution,  therefore,  requires  an  iterative  procedure.  The 
user  inputs  the  initial  guesses  for  the  constants  and 
iteration  then  continues  until  the  sum  of  the  partial 
derivatives  (Eqs  (A. 6),  (A. 7),  and  (A. 8),  is  less  than  10“6 
or  until  the  program  reaches  a  maximum  user  specified 
Iteration  count.  The  maximum  possible  iteration  count  and 
the  minimum  required  partial  derivative  sum  are  permanently 
set  at  150  and  10”6  respectively. 

The  foregoing  discussion  explains  the  procedure  for 
determining  the  coefficients  for  a  SINH  curve  corresponding 
to  a  particular  set  of  data.  Suppose  further  that  each  SINH 
representation  of  fatigue  crack  propagation  is  related  to 
other  SINH  curves,  the  relationship  depending  on  the 
differences  in  a  test  parameter  between  each  data  set. 
Consider  hold-time  as  the  changing  test  parameter  between 
each  data  set  as  an  example,  and  assume  that  the  points  of 
inflection  are  linearly  related,  namely 

^3fj  =  C33  +  C3n(C4,j)  (A. 13) 
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for  j  different  SINH  curves  which  are  fitted  to  j  different 
data  sets.  The  only  changing  parameter  between  each  set  is 
the  hold-time. 


Assume  also  that  C2  and  C4  are  related  to  hold-time  by 
Eqs  (A. 14)  and  (A. 15). 

c4,j  =  C35  +  CjgC-logCtfl))  (A. 14) 

c2,j  =  c37  +  C38(-log(tH))  (A. 15) 

The  functional  relationships  (-log(tH))  are  defined  in 
Table  1.  Coefficients  C33  through  C38  can  136  then 
determined  by  substituting  Eqs  (A. 13),  (A. 14),  and  (A. 15) 

into  Eq  (A. 9)  and  differentiating  with  respect  to  C33 
through  in  a  manner  analogous  to  that  used  in 

determining  C2,  C3,  and  C4  in  the  foregoing  discussion.  The 
constants  a  through  f  in  Eqs  (5)  are  the  appropriate  con¬ 
stants  C33  through  C38  of  Eqs  (A.13)»  (A. 14),  and  (A. 15). 
For  complete  development  of  the  SINH  model  see  (7).  For 
explicit  use  of  the  program  and  its  listing  refer  to  (16). 
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This  section  presents  the  computer  code  listing  of 
General  Electric’s  model.  This  listing  is  a  revised  version 
of  the  one  presented  in  (6).  The  program  is  written  in  a 
Fortran  computer  language  and  operates  from  the  Vax  11/780 
computer  hardware  system  located  at  the  Air  Force  Institute 
of  Technology,  Wright-Patter son  AFB,  OH.  This  program  is 
interactive  and  will  prompt  the  user  for  the  desired  test 
parameters.  After  calculation  of  the  six  coefficients,  the 
program  will  ask  the  user  if  values  of  crack  growth  rates 
are  desired.  If  so,  the  program  will  prompt  the  user  for 
the  required  information. 


program  ge  mod  i 

real  nl.n2.mm  | 

real  mrpt9,n 

open { 2 .file"’ ca lge2y * , status" ’ new’ ) 
open ( 3 ,f 1 le“ ’ ca 1 ge2x ’ , status" *  new ’ ) 

1100  print  *, ’Input  Temperature  (degrees  F):’ 
read  *,t 

print  *, ’Input  Stress  Ratio:’ 
read  * , r 

print  *, ’Input  Inverse  Freq.  (sec/cyc):’ 
read  *,tc 

print  *, 'Input  Hold  Time  (sec.):’ 

read  *,ht 

V025-0. 025-1 .e-20 

t999“ 1000 . - 1 .e-10 

tpt4  = .4-1 .e-10 

rpt9-. 9+1 .e-10 

s=-5.3+4.4e-3*t 

If (s. It. 1 .0)s=0.0 

n=8.e-04*t-0.04 

If (n. It. t .0 ) n= 1 . 0 

b=»  ( t- 1 200  >  *  .00 18-4.52 

lf(r.gt. . 5 ) xcons  = . 69897 

lf(r. le. .5)xcons2=.5 

1 f ( r . 1 e . .5>xcons=. 25527 

lf(r.gt. .5)xcons2=.  1 

If (b.gt.-4.52)b=-4.52 

yl=s*(alogl0(tc/tpt4) )**n+b 

n 1 = 1 . + 1 . 16e-02*(t-t999)-2.81e-05*(t-t999)**2 
sl=4.27e-13*(t-t999  )**4 .768+  .  2 
b 1 «7.607o-l 3* < i 400. -t> **4.605+4. 77 
bl"-bl 

y2*b 1  * ( a  1 og 10( tc / tpt4 ) )**nl+bl 

s2"2. I45e-04*(t-t999)**l .44 4 +.5 

b2»-2.628e-03*(t-t999)**l .051+6.201 

n2«. 9855+1 . 096e-03* ( t-t999 ) -4 . 366e-06* < t-t999 >«»2 
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b2=-b2 

y3«s2*(alogl0<tc/tpt4>  >**n2+b2 
s1ope=<yl-y2>/. 25527 

C=yl-slope*(alogl0(rpt9)-a1ogl0(l-r)> 

slopc2=(y2-y3)/.69897 

if<r.le..5)  go  to  10 

if(slope.lt.0.and.s1ope2.1t.0.and.3lope2.gt.slope.or.slope.gt.0. 

&a  nc).slope2.gt.0.and.3lope2.1t. si  ope) 

4c=y2-slope2*(a1ogl0( .5)-alogl0( 1-r)  ) 

10  continue 
dadn  = 10 . **c 
v  = 1 /tc 

1 f ( ht . eq . 0 . )go  to  590 
ef ct=0. 

tf(r.gt.-S)  go  to  530 

If (t. le. 1200. ) htb=0. 0065 *(t- 1000) 

efct=4.75e-03*< (alogl0( v/v025 )  )**. 9379-1.) 

If (t.ge. 1 200  )efct=efct*(t-l 200 ) 

If <t. le. 1200)ef ct=0 . 

If (t. le. 1 200 > hta =0. 005 *<t- 1000. > 

If (t.gt. 1 200) hta  =  1  . 

If (t.gt. 1 200 )htb~5.75e-04*(t-l 200 ) + 1 . 3 
htb=btb+ef ct 

htq=htb+hta*(a1oql0(ht/20) ) 

If (htq. lt.0. >htq=0. 
dadnl=10.**(htq+(alogl0(dadn) ) ) 
dadna=dadn 1 
530  continue 

etb=  1 . 456e-06* { t- 1000 . ) **2.275*1  ( a  1 og 10( v/ v025 > )**<9. 15e-04 
4* 1 1- 1000 . ) ♦ . 652  )  - 1  .  ) 
br 5=3 . 809e-02* { t- 1000 . ) ** . 5905+etb 

evt=( .55-1 .739e-06*( t- 1000.  >**1 .998>*( ( a  1 og 10< v/ v025 ) )**.25-l . ) 
cr5=.55+7. 1578e-03*( t- 1000 . )-l . 55e-05* ( t- 1000 . >**2.+evt 
htq5=br5+cr5*(alogl0{ht/20. ) ) 

I f ( htq5 . 1 t .0 ) htq 5=0. 
dadn5=10.**(htq5+<alogl0(dadn> ) ) 

If ( r . le. .5  )dadnb  =  dadn5 
If ( r . gt . .5)dadna=dadn5 
lf(r.le..5)  go  to  535 

evt9=<0.G4-3. 167e-05*(t- 1000 . >**1 . 5 1 3 ) * ( ( a  1 og 10( v/ V025 ) ) ** . 32- 1 . > 

c259=.64  +  6.572e-03*(t-l 000 ) -  1 . 4  8e-05* ( t- 1000 ) **2 . 

cr9=c259+evt9 

br9=. 1 5+. 3543* (t- 1000. >**.1339 

xz  =  ( . 13  +  4. 9e-05*(t- 1000. )**1 .5G>*< ( a  1 og 10( v/ v025 ) >**1.8-1.) 

br 9=br9+xz 

if (cr9. It.0)cr9=0. 

htq9=br9+cr9*(alogl0( ht/20. ) ) 

If (htq9. lt.0. ) htq9=0 . 
dadn9=10.**(htq9+lalogl0(dadn> ) ) 
dadnb=dadn9 

535  dadn=(alogl0{dadnb)+{<alogl0(dadna/dadnb)>/xcons)*(alogl0((l-r) 
4/xcons2  )  )  ) 
dadn=10.**dadn 

590  dl=l . 25618-2. 83Se-09*l t-t999)**3. 10561 
d2-l .37183-5.G3102e-10*(t-t999)**3. 30553 
a=2. 081-1 dl*<a log 10(rpt9/( 1-r ) ) )**d2> 
mm« . 1 3 

If (r . le.0.5)go  to  17 
mrpt=9.898e-06*( t-t 9991**1. 5847 
s 1 pm= ( . 13-mrpt9)/. 69897 
mm* . 13-s1pm*(alogl0( .5/( 1-r ) > ) 

17  continue 

de 1 k  =  10 . **a*dadn  **mm 
xkc= 1 22* (  1-r  ) 

xks-10.**(  1  .-{  (  .  232  +  8. 406e-03*<t- 1000)**.  729  )*(a  log  101  .9/U-r  >  )  )  1  ) 
1f(delk.gt.xkc)delk=xkc- .0001 
vht= 1 / ( v  +  ht  > 

I f ( r . gt . . 5 )go  to  600 
ttemp=t 

I f  <  t . gt . 1 300) t temp  =  1 300 

e= 1 . 55e-02*  <  ttemp- 1 200 ) * ( ! a  1 og 10( v/ v025 ) >**.0902-1 ) 

If (e. lt.0. ) e=0 . 

pr im=2. 7+6. 88e-03*(t- 1000)- 1 . 94e-05* < t- 1000 > **2 . +e 
If (ht.eq.0. )go  to  599 


a25t=2 . 737e-07* 1 ttemp- 1000) **2 . 7  4 
a vt=a25t+e 

bvt=5 .34Ge-03*l ttemp- 1000)** .7331 
qvtht=avt+bvt*(alogl01ht/20. >  > 
pr im=pr im-qvtht 

if (t. le. 1200) up =2. 7-1 . e-03* < t- 1 000 > 
if It.gt. 1200) up =2. 6-. 13 72* It- 1200) **.5 194 
if(pr im. lb. up )pr im=up 

599  pr  i  ma  =pr  i  in 

600  continue 

gama  =  1 . 7-8 . 25e- 10* < t- 1000 ) ** 3 . 54 

e= 1 .734e-10*<t- 1000 ) **3 . 807* ( (alogl0! v/v0 25) >**gama-l ) 

pr im=2.25+6.875e-03*(t-1000>-2. 1 88e-05* ( t- 1000 ) **2 . +e 

if  1 ht .eq . 0)go  to  610 

a25t  =  1 .652e-02*(t  -  1000 )  **  .  66  1  6 

alpha*. 5223  +  4. 23e-03*< 1 400- 1  > 

iflalpha.gt. I ,37)alpha=l .37 

beta*. 152+7. 165e-03*< 1400-t) 

iflbeta.gt.l .59)beta=l .59 

b25t= . 1333*(t-l 000 ) *  * . 1 1 1 

a vt=a25t+3 . 1 3e-03* 1 t- 1000 ) **2 . SG* 1 < a  1 og 10( v/ V02G > )**alpha-l ) 
bvt  =  b25t  +  4 . 53e- 1 i *  <  t- 1000 ) **3 .755*1 (alog!0( v/v025  > )**beta-l ) 
qvtht=avt+bvt*la1ogl0(ht/20)) 
pr im  =  pr im-qvtht 

up  =  2.25-2.613e-09*(t- 1000  >**3.426 
610  i  f < pr 1 m . 1 1 . up ) pr i m  =  up 
iflr.le..5)primb=prim 
i  f  1  r  .  g  t .  .  5  )  p  r  i  m  a  *  p  r  i  m 
i f 1 r . 1 e . . 5 ) go  to  620 

e  =  4. 351 e-03* It- 1000) **.7225*1 ( a  1 og 101 v / V025 ) >**1.34-1) 
pr im=l .4+5.425e-03*(t- 1000  > - 1 . 963e-05* I t- 1000 > * *2 . +e 
I  f ( ht . eq . 0 . ) go  to  619 
a lpha=l .705+2. 66e-03*( 1400-t > 
a25t*7.87e-06*It- 1000 ) ** 1 .716 
if (alpha.gt.2.237)alpha=2.237 

avt=a25t+7.S71e-0G*l  t-1000)**l .716*1 1  a  1 og 101 v/v025 > )**alpha-l ) 
b25t= 1 .575e-03*(t - 1000 )-2.875e-06*l t- 1000 ) **2 . 
gama=. 81 94-2. 49e-03*< 1400-t) 

I  f 1 gama . 1 1 .0. 322 ) gama= . 322  i 

evt=b25t+ 1 Ialogl0(v/v025) )**gama-l )  1 

b vt=b25t+evt 

qvtht=avt+bvt*la1ogl0(ht/20) ) 

pr  Im=pr  im-qvtht 

up* 1 .5-7. e-03* It- 1200) 

iflt.lt. 1 200 ) up= 1 .4+1 .e-04*lt-1000) 

If  1 pr im . 1 t . up )pr im=up 

619  pr i mb*pr 1 m 

620  continue 

prlm=lalogl0{primb)+i(a1ogl0(prima/primb!)/xcons)*(alogl01il-r> 
Jr/xcons2 ) ) ) 
dadnp* 10. **pr i m 

ape=a log  101 2. 5) +2. 41 e-03* It- 1000) 
q-10.**Iape+.57*lalogl01 . 1/r ) ) ) 
if  1 q . gt . 3 .  )q*3 . 

srd=lsqrtlq)*Ialog!xkc)-alogldelk)))/lalogldelk)-aloglxks)) 
d*- I srd**2 . ) 

p=dadnp-lq/lalog(delk)-aloglxks) ) )+Id/IalogIxkc)-alog(delk) ) ) 
bp*a log! dadn )-q*a log! a  log! del k ) -a  log! xks ) ) -d*alog I  a  1 ogl xkc ) 
c-a 1 og 1  del k ) ) 

b=bp-p* ialog(delk)-aloglxks)) 
print  40,b,p,q,d,xks,xkc,dadn,dadnp 
print  4 1 . t , r , v , ht 

40  format(6(2x,f8.4) ,el4.5,2x,f8.4> 

41  f ormat 1 4 1 f 10. 4 . 2x ) ) 

print  *, ’Desire  da/dn  values?  1  for  yes,  0  for  no' 
read  * , f lag  1 

i  f 1 f  1  agl .eq .0 )  go  to  999 

1000  print  *, 'Input  desired  delta  K  Iksi  root  In)’, 
read  *,dk 

dadno  =  exp 1 b  >  *  1  I dk/xk s ) **p ) *  1 1  a  1 og 1 dk /xk s ) ) **q  )  * 

4 1 1  a  1 og 1 xkc/dk ) ) * *d ) 


.  -  m  _  Vr**  m  «  . 


Afifiendii  C 


AF115  is  a  gamma-prime  strengthed  nickel-base  alloy  in 
which  titanium,  aluminum,  columbium,  and  hafnium  are  gamma 
prime  formers  and  chromium,  cobalt,  molybdenum,  and  tungsten 
are  strengtheners  of  the  gamma  matrix.  The  AF115  powder  was 
produced  by  vacuum  melting  from  virgin  material  and  argon 
spray  atomized  to  powder.  Results  of  the  chemical  analysis 
conducted  by  the  General  Electric  Company  (6:2)  of  the 
powder  are  shown  in  Table  4. 

This  composition  was  produced  by  Carpenter 
Technology  (6:7).  They  produced  a  cylinderical  log  with  the 
dimensions  shown  in  Figure  22a  along  with  the  location  of 
the  14  slices  used  to  produce  test  specimens.  Figure  22 b 
illustrates  how  the  compact  tension  specimens  were  machined 
from  each  slice.  Seven  of  the  slices  (1,  2,  3,  6,  12,  13, 
and  14)  were  used  for  density  characterization,  stress 
qualification,  and  thickness  evaluation  (6:3-14).  The 
remainder  of  the  slices  each  produced  8  compact  tension 
specimens  which  were  used  for  cyclic  testing.  The 
configuration  of  the  machined  compact  tension  specimens  is 
illustrated  in  Figure  23.  Because  of  the  undersize  diameter 
of  the  log,  these  specimens  contained  one  beveled  corner. 
It  was  located  near  one  of  the  loading  holes  and  had  no 
influence  on  crack  growth  results  (6:6). 


The  compact  tension  specimens  were  precracked  at  either 
room  or  elevated  temperature.  Initially,  the  specimens  were 
precracked  at  room  temperature  at  30Hz  with  cyclic  loads 
being  stepped  down  in  10$  increments  as  a  function  of 
precracked  length.  The  final  precrack  length  and  applied 
loads  produced  a  stress  intensity  range  (aK)  less  than  or 
equal  to  the  initial  stress  intensity  of  the  subsequent 
elevated  temperature  cyclic  crack  growth  tests.  About  half 
of  the  specimens  were  precracked  in  this  fashion.  During 
the  early  stage  of  testing  it  appeared  that  the  precracking 
was  interfering  with  the  normal  crack  growth  behavior.  The 
curvature  of  the  crack  would  rapidly  change  shortly  after 
initiation  of  testing.  Because  of  this  effect,  precracking 
was  conducted  at  elevated  temperatures  within  the  system 
used  for  testing  so  that  the  final  precracking  was  done  at 
the  test  conditions  loads.  Even  then  a  few  of  the  specimens 
produced  uneven  crack  fronts  (6:18). 

Specimen  heating  was  provided  by  a  specially  designed 
split  shell  three  zone  resistance  furnace  equipped  with 
quartz  viewing  windows.  Conventional  20X  traveling 
microscopes  were  used  to  monitor  crack  growth  on  both 
surfaces  of  the  specimens. 

At  the  conclusion  of  each  test,  the  raw  crack  length 
versus  accumulative  cycle  data  were  reduced  to  cyclic  crack 
growth  rate  (da/dN)  by  use  of  the  seven-point  sliding 
polynomial  technique  recommended  by  the  American  Society  for 
Testing  and  Materials  (ASTM)  (17).  The  cracks  lengths  were 


adjusted  for  curvature  by  the  ASTM  recommended  quarter  width 
average  crack  length  and  maximum  crack  length 
technique  (18).  The  stress  intensity  range  value  associated 
with  the  mid-crack  length  of  each  set  of  consecutive  seven 
points  was  determined  by: 


AP  2  +  a/W 

ak  = - —  - —  (0.866  +  4.64(a/W) 

B(W)1/2  (1-a/W)3/2 


-  13.32(a/W)2  -  5.62(a/W)3) 


(C.1) 


where 


aP  =  change  in  load  (Pmax  -  Pmin) 
B  =  specimen  thickness 
W  =  specimen  width 
a  =  average  crack  length. 


Thus  pairs  of  da/dN  and  aK  values  were  available  for  model 
evaluation. 
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